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ABSTRACT OF THE DISSERTATION 


A new classification method for arrhythmia detection for applications in the automotive industry 
By 
Natalia Estefanía Jurado Espín 
Automotive Engineering Masters in MIR 09 
ITESM, Toluca, 2015 


Dr. Luciano Chirinos Gamboa 
Dr. Jaime Salazar 
Dr. Alfredo Santana 


In this study, the design, validation and implementation of an artificial neural network (ANN) solution for 
arrhythmia detection is presented. Four types of ECG patterns were chosen from the MIT-BIH database 
to train the network, posteriorly tested with real time monitored heart beat electrical signals. Compared 
with other works, this research significantly reduces complexity without compromising its overall accuracy 
by feature reduction techniques. The ANN was implemented on Android"M based mobile devices (e.g. 
smartphones and tablets) and includes a QRS window extractor, a classifier and further blocks to detect 
abnormal heart beats with an overall accuracy of 91.7%. Additionally, the system generates messages to 


alert the user and helping to develop some future automotive safety applications. 


1x 


Chapter 1 


Introduction 


Heart abnormalities are one of the main road death causes, moreover, sudden cardiac death as a result of 
an arrhythmia is the most feared complication among drivers. To monitor these conditions and to prevent 
an accident must be a priority for the security of the driver. In this sense,the use of on-board electronics 
and in-vehicle information systems give us the opportunity to collect and analyze physiological data during 
driving tasks intended to help prevent casualties by stopping the car if vital life signs of the driver are not 
well. For this purpose, the use of the computational power of daily artifacts such as smartphones or tablets 
can be used for monitoring and linking the patient with the car, representing an implementation challenge 


that is object of the study in this thesis. 


1.1 Research Scenario 


Road traffic deaths caused around 1.25 million of deaths each year [24] (plateaued since 2007), from this 
total 167,140 happened in Mexico (Figure 1.1), being the 12th most affected country in the world. Children, 
pedestrians, cyclists and older people are among the most vulnerable of road users (Figure 1.2). In the face 
of rapidly increasing of this problem, heads of world are taking actions. For instance, in September 2015, the 
United Nations General Assembly adopted the historic Sustainable Development Goals (SGDs), where one 
of them is to halve the global number of deaths and injuries from road traffic crashes by 2020. Moreover, 


governments are initiating programs of medical examinations for all drivers. 
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Figure 1.1: Road Traffic Deaths by Country 


The roots of road accidents are being studied and there is also foundations to say that chronic medical con- 
ditions and traffic accidents is a subject of controversy. From one side, an estimated of 17,32 million people 
die every year due to cardiovascular diseases, being the top cause of premature death: Coronary (ischaemic) 
heart disease. On the other side, according to the Driving and Heart Disease Journalist [5] sudden cardiac 


death as a result of an arrhytmia is the most feared complication among drivers. 


Over the last year, around 100,000 people die and 50 million people are hurt [24] died in a road 
accident caused by a heart condition. This number represents the 7,1% of the total deaths while driving 


1.1. To reduce road traffic deaths is clearly necessary to work in partnership: governments, local authorities, 


Death Cause Percentage 
Epilepsy 38% 
Blackouts 21% 


Diabetes and insulin | 18% 
Heart Condition 8% 
Stroke 7% 
Others 7% 


Table 1.1: Road death causes, World Health Organization 2015 [24] 


police forces, individual road users and, above all, car manufacturers. Safe vehicles play a critical role in 


averting crashes and reducing the likelihood of serious injury. Over the past few years a combination of reg- 
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Figure 1.2: Road Accidents by Road User 


ulatory requirements and consumer demand has led to increasingly safe cars. On this subject, automotive 
companies are to join efforts to achieve effective and long-lasting improvements in road safety. In order to 
decrease car accidents due to heart conditions, electrocardiogram (ECG) signal acquisition and analysis has 


been implemented in vehicles; but this systems had not been fully commercialized. 


1.1.1 Previous Works 


Irregular heartbeat (arrhythmia) always precedes a heart attack and it is an accurate indicator of any heart 
condition. Arrhythmia can be detected by scanning the ECG signal, analyzing it, making decisions in order 
to control the car and even sending patient information by car connection with hospitals and ambulances. 
Diagram 1.3 shows the whole intended system explained by modules and subsystems. Such examples have 
been object of study in different research and innovation centers in the automotive industry, the most relevant 


examples are: 
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Figure 1.3: Arrhythmia Detector System 


Ford Safety Seat 


A heart attack behind the wheel can render much of your car’s safety equipment moot despite decades of 
advancement. Ford Motor Company has tasked its European Research Center in Aachen, Germany with 
finding a way to reduce accidents caused by drivers experiencing heart trouble. According to the automaker, 
their prototype seat with contactless electrocardiogram technology can warn drivers to seek medical attention 
immediately by scanning for potential cardiovascular trouble through clothing. 

According to Ford, by 2025, nearly one quarter of Europeans will be at least 65 years old, a figure that 
increases to one third by 2050 [9]. As risk of cardiovascular disease increases with age, drivers suffering from 
conditions like angina can be more than fifty percent more likely to be involved in an accident, so there’s a 
lot of motivation to push this technology. 

The prototype seat presently delivers accurate readings for 95 percent of drivers 98 percent of the time. 
Ford researchers are now looking into how the seat can be integrated with the other safety systems to work 


together to protect drivers who experience issues behind the wheel. 


Toyota ECG Sensors in Steering Wheels 


Toyota is working on a multitude of technologies that will make driving safer. A pre-collision system with 


a steering-control feature, a pop-up hood and blocking of high-beam headlights, and also an ECG detecting 


Controls 


Figure 1.4: Ford plus Medtronics for In-Car Health Managment 


driving wheel system has been studied in their facilities in Japan. 

Whereas Ford is opting for a contactless ECG-sensing driving seat, Toyota is going the more obvious way, 
by direct contact sensors in the steering wheel. Of course this depends on the driver keeping hands on the 
steering wheel long enough to detect any abnormal heart rhythm, which might not be the case in more acute 
conditions such as VF. An optical sensor mounted in the steering wheel picks up a single-lead ECG signal, 
and alerts the driver through the display about any abnormalities that need attention. In an experimental 
system installed in a Prius [32], the chart of the drivers heartbeat and electrocardiogram can be viewed on 
the display of the in-car navigation system. In addition to alerting of any serious conditions, Toyota also 
sees the system being used as a sort of daily health check. 
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Figure 1.5: Toyota to Integrate ECG Sensors Into Steering Wheels 


BMW Emergency Stop Assistant 


The BMW Emergency Stop Assistant is based on two types of sensors. On one hand it has sensors designed 
to monitor the vital life signs of the driver and on the other hand it benefits from BMW’s ConnectedDrive 
assistance systems [2]. When the vital signs monitor detects the driver has become impaired, the Assistant 
kicks in and automatically activates an autonomous driving mode. The system then turns on the hazard 
warning lights and starts maneuvering, while taking into account the traffic conditions in which it finds itself 
in. When it detects it has reached the outer edge of the road, the vehicle grinds to a halt. Once the safety 
of the car and of the driver have been insured, an emergency signal is sent to rescue services, containing the 


data required for emergency services and traffic coordinators (should traffic be disturbed by the incident) to 


move in. 


In a nutshell, these three works show the evolution of On-Road ECG Tech with one particular target: To 


Figure 1.6: BMW Automated Automobile 


reduce the number of road deaths caused by heart conditions. 1.2 


Company 


Challenge 


Achievements 


Next Step 


Ford 


Contact-less ECG technology 


98% of the time over 95% 
of the drivers of accurate 
readings 


Incorporate a User Interface 


Toyota 


Incorporation of a User Interface (UI) 


Contact-less ECG sensors 
installed in the steering 
wheel connected to a nav- 
igation system 


Road assistance based on ECG measurements 


BMW 


On Road Assistance 


Vital Life Signs plus Con- 
nectivity plus On Road 
Assistance 


Implementation and commercialization 


Table 1.2: Previous Works. 


1.2 Motivation 


In order to reduce the risk of accidents, detecting arrhythmia on time can buy some vital minutes to the 
driver and any other road users. Most of the previous works and efforts in the industry show the need of 
specific computers in order to accomplish the task of analyzing ECG signals. The implementation of the 
system using specific modules of hardware and software increase the cost significantly. The challenge is to 
find a cheap and trustful implementation which can be a realistic choice not only for the driver but also for 


the car manufacturers. 


First, there were hands-free devices, then voice-directed gadgets. Now the auto industry is finding dif- 
ferent solutions for particular type of road users taking advantage of the widely use of smartphones and 
tablet.Technology is certainly moving in that direction and also offers great chance to reduce costs by mak- 


ing use of the computational power of devices previously acquired. 


Jen-Hsun Huang, the co-founder and executive director (CEO) of the NVIDIA, declared: In the future, were 
going to have all cars be connected, and they will be connected to the cloud. You will get software upgrades 


over time and, as a result, your car gets better and better over time. 


Moving to a software defined car is happening right now with the help of all the components of the automotive 
industry world. Talking about connectivity, one can dream of self-drive cars that can make vital decisions 
when preventing road trafic deaths. All, cars, ambulances, family doctor and patient (driving cars with no 
risk) can be connected constantly ready for any possible inconvenience. Therefore, this study is motivated 


by the facts that: 


e Heart condition is a major problem for car drivers. 


e Early detection of heart conditions (e.g. arrhythmia) through ECG sensing can help to prevent fatal 


car accidents. 


e Mobile data networks availability can help us to communicate with systems and people (almost in real 


time) needed to assist the driver. 


e Most of previous works are not suitable for implementation on mobile devices, new algorithms need to 


be developed to address the limited processing capabilities. 


e Addressing solutions that have to be implemented with dedicated hardware increases the cost signifi- 


cantly. 


e Implementing the software taking into advantage the existing hardware reduces costs and also let the 


drivers choose the modules more suitable for their health profile. 


There is one major obstacle in achieving a reliable system in terms of reading one's vital signs. One has 
to find a way for the sensors to detect the physiological data of the driver (heart rate, breathing) without 
entangling the drivers in wires. Another issue is the development of an algorithm which will allow emergency 
services to pinpoint the vehicle within its lane, as well as make the car’s computer able to identify objects 


in the area surrounding the car. 


As separated in subsystems showed in Figure 1.3, each challenge has to be studied separately. But finally, 


the intended system is showed in Figure 


1.3 Problem Statement 


In this study I intend to find out if it is possible to implement an arrhythmia detection system, for preventing 
car accidents, on a mobile platform. To achieve this target, the challenge is to develop a new classification 
method for arrhythmia detection suitable for cheap and trustful implementation for future automotive ap- 
plications. In addition, the system must be easy to implement in terms of algorithm complexity, possible to 
connect with different car user interfaces and reliable to make vital decisions. 

The boundaries of this proposal are that the classifier will be based on a Back Propagation Training method 
trained with standardized data from MIT-DB, the three most common types of arrhythmia tachycardia, 
atrial flutter and fibrillation will be taken into account for this study and tested and implemented on a 


mobile platform (Android). 


y 
v 


Figure 1.7: Intended System 


1.3.1 Objectives 


General 


To create an algorithm for classification of arrhythmia based on a new set of QRS complex features for 


efficient implementation on a mobile device. 


Specific 
e To digitally process the ECG signal for proper detection of the QRS complex vector. 
e To train a Back Propagation Neural Network to classify different arrhythmia types using MATLAB. 


e To explore feature classification alternatives to reduce the complexity of the implementation without 


sacrificing the accuracy. 
e To test the trained network with new data in online application. 
e To implement a feed-forward network for arrhythmia classification on a mobile platform. 


e To develop some experiments to show the performance of the algorithm 


1.4 Contributions 


Making cars safer is a critical component of saving lives on the roads, in this sense, this study granted two 


main contributions: 


e A 61,9% of the feature complexity was reduced compared to A Recognition of ECG Patterns Using 
Artificial Neural Network study, which was a referent for this current research. In order to find the 
best combination of features, a series of repetition was selected based on the extraction of main char- 


acteristics of the ECG signals and QRS complexes. 


e An Android application was implemented on a tablet and smartphone in order to classify the samples 


in a mobile platform. 


e In the Automotive Mechatronics Research Center, this study is a pioneer work and leads a possible 


future series of investigation of improving driver’s safety on road. 


1.5 Content Description 


After getting deep into the research scenario and previous works. The motivation which justifies the study 
of this topic is also stated in this section. Each specific target conducting to the general objective are pin- 
pointed to state the contributions of this research. Chapter 2 contents the background of this study. The 
arrhythmia concepts and nature, ECG measurements and data bases, Neural Networks (NN) fundamentals 


and the mobile platform's operating systems. All of the necessary concepts required to understand the 


10 


methodology, results and conclusions as a contribution of this study. Chapter 3 describes the methods, step 
by step, followed over this work in order to select the combination of features, train the neural network and 
create the algorithm to be implemented in a mobile platform. Chapter 4 states the results obtained as the 
specific objectives lead to the main target.Chapter 5 concludes the work and also states the future work in 


automotive applications. 
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Chapter 2 


Background 


From the ECG signal acquisition, QRS complex characteristics, Artificial Neural Networks (ANN) and 
Android Mobile Platform, all the necessary knowledge for understanding the intend of the present work is 


described. 


2.1 Generation of ECG signal 


2.1.1 Anatomy of the heart 


The heart consists of four major regions: left atrium, right atrium, left ventricle and right ventricle. Left 
and right parts of the heart are separated by the septum [30]. In between atria and ventricles there are 
valves that enable blood flow from the atria to the ventricles but inhibit backflow to the atria. In between 
ventricles and arteries there are valves, too, that inhibit blood backflow from the arteries to the ventricles 


after the contraction of the heart. 


2.1.2 Pumping Function 


The right ventricle is responsible for pumping blood to the lungs, where the blood is oxygenated. From 


the lungs the blood returns to the heart and it flows into the left atrium. With each heart beat, the atria 
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Figure 2.1: Intended system in terms of signals as a system 
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Figure 2.2: Anatomy of the heart 
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contract and discharge their content into the ventricles. Shortly after the atria's contraction, the ventricles 
contract, too. Blood from the left ventricle is pressed into the aorta and from there it is disposed to the 
various tissues of the body. From the body the blood is reunited within the vena cava superior and the 
vena cava inferior, that lead to the right atrium, which contracts prior to each ventricular contraction and 
discharges its blood to the right ventricle. The contraction phase of the heart is called systole, and it is 
followed by a period of relaxation, where the atria and ventricles are refilled with blood [14]. This phase is 


called diastole. All of this process can be seen in 2.3. 


2.1.3 ECG Signal 


An electrocardiogram signal (ECG) contains a large amount of information which can be used for analyzing 
many attributes of the electrical activities of the heart. Typically, an ECG varies from person to person 
due to the difference in the position, size, anatomy of the heart, age, body weight, chest configuration, and 
various other factors. An ECG is a representative signal of cardiac physiology, which can be very useful in 
diagnosing cardiac disorders. The ECG was originally observed by Waller in 1899 [30], [18] and [14]. In 
1903; Einthoven introduced the electrophysiological concepts still used today, which includes the labels of 
the waves. To this day, an ECG is the main way of gathering information from the body in order to better 
analyze the hearts activities 2.3. The heart is one of the most important organs in the body due to the fact 
that it supplies the body with oxygen. 

The corresponding part in the heart plays their respective roles. Sinoatrial (SA) node will excite the beats 
that caused heart muscles to contract. Below is shown the location of Atrioventricular (AV) node and 
SA node which are responsible for generating ECG signal in the humans heart. Figure 1 shows the basic 
structure of heart. The contraction of hearts muscles soon will be recorded as an electrical activity of the 
heart called ECG signal.Based on the pattern of ECG recording, heart status could be identified whether 
possessed of any cardiac arrhythmias or otherwise. As known, heart muscle possessed the characteristic of 
depolarization and repolarization. Depolarization is referred to the electrical potential activity excited by 
heart muscles while repolarization is a relaxation state when the heart changing back to its original position. 
P wave generated due to the atrial depolarization, QRS complex represented ventricular depolarization while 
T wave represented ventricular repolarization. Figure 2.4 shows the corresponding part of heart function 


with respect to the ECG signal obtained. 
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Figure 2.3: How ECG is generated from the heart 


Figure 2.4: Diagram of ECG signal outputted by the heart 
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2.2 Heart (ECG) abnormalities 


Abnormalities happened in the respective waves segment will provide ideas to doctors and cardiologists at 


where the part of heart is having problems. 


Types of Arrhythmia 


1. SVT (supraventricular tachycardia) is a heart arrhythmia that is due to a fast heartbeat in the atria. 
This can occur in patients for a few minutes to a few days and can cause dizziness, rapid breathing, 
and chest pain. SVT is rarely life threatening but it can be treated by a few medications or physical 


maneuvers that are relatively simple. 


2. VT is a heart arrhythmia that originates in the ventricle of the heart. Unlike SVT, this is a life 
threatening arrhythmia that can cause ventricular fibrillation and sudden death. 80 % of cases of 


sudden cardiac death (SCD) are caused by spontaneous VT. 


3. BBB is due to disruption in the normal flow of electrical pulses. Typically, BBB has QRS width 
duration greater than 120 ms. Right and left BBB are possible and the ECG signal will be shifted left 


or right towards the P or T waves respectively. 


4. Bradycardia is defined as a slow heartbeat below 50 beats per minute (BPM). This can be caused due 


to drug use or heart disease. 


5. Tachycardia is defined as a fast resting heart rate of above 100 beats per minute. Like Bradycardia, 


this is quite similar other than the heart is racing due to bad diet or poor heart condition. 


6. Premature atrial contractions (PACs), also known as atrial premature complexes (APC) or atrial 
premature beats (APB), Are a common cardiac arrhythmia characterized by premature heartbeats 
originating in the atria. While the sinoatrial node typically regulates the heartbeat during normal 
sinus rhythm, PACs occur when another region of the atria depolarizes before the sinoatrial node and 


thus triggers a premature heartbeat. 


7. Accessory pathway tachycardias (bypass tract tachycardias) are a rapid heart rhythm due to an extra 
abnormal pathway or connection between the atria and the ventricles. The impulses travel through the 
extra pathways (short cuts) as well as the normal AVHIS Purkinje system. This allows the impulses 


to travel around the heart very quickly, causing the heart to beat unusually fast. 
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Figure 2.5: The most common arrhythmia 


8. Atrial fibrillation, a very common irregular heart rhythm. Many impulses begin and spread through 
the atria competing for a chance to travel through the AV node. The resulting rhythm is disorganized, 
rapid, and irregular. Because the impulses are traveling through the atria in a disorderly fashion, it 


results in loss of coordinated atrial contraction. 


9. Atrial flutter, An atrial arrhythmia caused by one or more rapid circuits in the atrium. Atrial flutter 


is usually more organized and regular than atrial fibrillation. 


10. Premature ventricular contractions (PVCs) are early extra beats beginning in the lower chambers of 
the heart (ventricles). PVCs are common.Most of the time they cause no symptoms and require no 


treatment. In some people, PVCs can be related to stress, too much caffeine or nicotine, or exercise. 


Most common arrhythmia 


Atrial fibrillation is the most common sustained arrhythmia, increases with age, and presents with a wide 
spectrum of symptoms and severity. Paroxysmal, persistent, and permanent forms require very individualized 
approaches to management [33]. The most common arrhythmia presented by this referenced document are 


shown in figure 2.5 
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Figure 2.7: Synapses 


2.3 Neural Networks 


A neural network is a type of computational model which is able to solve multi problems in various fields. It 
processes the information in a similar way as the human brain concept processing the information. Basically, 
neural network consists of large processing elements called neurons working together to perform specific 
tasks. As in the human brain, there are thousands of dendrites which contain information signals shown in 
figure 2.6. They transmitted the signals to the axon in the form of electrical spikes. The axon then sends the 
signals to another dendrites causing to a synapse shown in figure 2.7.This synapse occurred when excitatory 
input is sufficiently large than the inhibitory input, and this concept of signal transmission also depicted 
on how neural network process inputs received [16]. During processing stage, inputs will be trained to 
produce desired target outputs until it come to a generalization stage. Generalization means at a condition 
where the network is able to recognize the inputs and the corresponding targets after undergone the training 
given. Then, the network will be tested by given new inputs signal to evaluate its performance and to see 
how accurate the output produced will be comparing to the target. Figure 2.8 and 6 depicted the analogy 


of human brain concept process the information to the neural network system [25]. Neural network consists 
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Figure 2.9: Neural Network 


of several architectures, from simple structure until the complicated ones. 


2.3.1 Single-Layer Feed forward Networks 


This is the simplest form of network architecture with only single layer output without any hidden layer. 
An input layer of source nodes will directly projects onto the output layer of neurons or computation nodes, 
just in one way but not vice versa. Single layer is referring to the output layer which is just single output 
and not considered the input layer of source nodes since no computation performed there. Figure 2.10 is 


shown the corresponding network. 


2.3.2 Multi-Layer Feed forward Networks 
This second layer is differing from above since it has one or more hidden layers. The computation also takes 


place in these hidden nodes. The hidden nodes are also used to intervene between the external input and 


the network output with respect to the networks manner. Figure 2.11 shows the network architecture. 
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Figure 2.11: Multi-Layer Feed forward Networks 
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Figure 2.13: Recurrent network wit feedback 


2.3.3 Recurrent Network With No-Self Feedback Loops and No Hidden Neu- 


rons 


A recurrent neural network has at least one feedback loop. It may consist of a single layer of neurons 
with each neuron feeding its output back to all the input neurons as illustrated below. The feedback loops 
increased the learning capability of the network and on its performance. Besides, these feedback loops are 
also associated with unit delay elements (z-1) which result in a nonlinear dynamical behavior in a condition 


when neural network contains nonlinear units. Figure 2.12 shows recurrent network without feedback [1]. 


2.3.4 Recurrent Network with Hidden Neurons 


This structure distinguishes itself from subsection 2.3.3 with hidden neurons. The feedback connections 
originate from the hidden neurons and also from the output neurons. The structure illustrated in figure 
2.13. In neural networks, the input layer is passive while the hidden nodes and output layer are active and 


normally been activated by transfer function. These active nodes will modify weight and bias values to 
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an optimum number where the network is best works at. There are many activation functions which can 
be applied such as radial basis (radbas), competitive transfer function (compet), positive linear (poslin), 
saturating linear transfer function (satlin) and many more. The most commonly used is including hard- 
limit transfer function, linear transfer function and log-sigmoid transfer function [30],[14] and [10]. Neural 
networks need to be trained with suitable learning algorithm training functions corresponding to a network 
type. Among of the training functions are gradient descent backpropagation (traingd), Levenberg-Marquadt 
backpropagation (trainlm), gradient descent with adaptive learning rule backpropagation (traingda), random 


order incremental training with learning functions (trainr). 


2.4 Review of Neural Networks implemented for Detection of Car- 


diac Arrhythmia 


The result of different neural networks cannot achieve 100 % or very high percentage of accuracy since small 
data set was used. Thus, some recommendations could be suggested to improve networks performance in 
order to obtain approximately accurate outputs. Another suggestion, with high complexity implementation, 
must be to increase the numbers of hidden neurons and retrain the network several times. Furthermore, the 
usage of larger data set can be tried, so that the network will learn more and expose to enough training. To 
try a different training algorithm as well as adjust the initial weights and biases to new values. Then, train 
the networks again for several times until it reaches desired target. These improvements are recommended 


to be done in the future for better result by [21] and taking into account over this study. Table 2.1 shows a 


comparison of different neural networks with different algorithms . 
Neural Network Data Analized Accuracy | Approach 
SOM 500 Data Samples 85 % | FFT 
BP 49 data samples 81.6 % | Pattern Recognition for feature extraction 
Fuzzy Kohonen Network 400 data samples 86.21 % | Multifractionality 
BP Near about 60 samples 96 % | Empirical mode decomposition (EMD) 


Table 2.1: Comparison of most common Neural Networks used for arrhythmia detection 
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2.5 Mobile Implementation for Arrhtythmia Detection Efforts 


Mobile devices like smartphones and tablet computers continuously grow in processing power and become an 
integral part of daily life, even in developing countries [29]. Recently, such mobile devices are also used for 
biomedical signal processing and ECG analysis [12], [7], [8]. In 2011, PhysioNet/Computing in Cardiology 
[29], [13] arranged a challenge to develop efficient algorithms to improve the quality of ECG recordings using 
mobile devices and to improve efficiency of ECG diagnosis. Scully et al. [28] demonstrated that physiological 
features like heart rate, breathing rate and blood-oxygen saturation can be extracted using camera recordings 
from mobile phones. The authors concluded that all processing could be performed on modern mobile 
devices. However, they did not share their software implementation. Furthermore, algorithms for arrhythmia 
detection were developed and assessed that were capable of handling the hardware restrictions of mobile 
devices. However, they needed medical expert intervention in order to operate. In summary, the previous 
work in the literature focused on developing algorithms for ECG monitoring either did not share an actual 
software implementation, did not specifically evaluate arrhythmia detection on mobile devices, or needed 
expert intervention in order to operate. Our contribution in this paper is threefold. Firstly, we provide an 
algorithm for real-time detection of QRS complexes and automated, intervention free normal /abnormal heart 
beat classification, which extends well-known analysis methods [26]. Secondly, we present an implementation 
of the algorithm in an Android-based ECG monitoring applications, which can process ECG signals in real- 
time by accessing a ShimmerTM sensor node via Bluetooth or by using a database of pre-recorded data. 
Thirdly, we conduct a detailed evaluation of the application and the implemented algorithms in respect to 
QRS detection and abnormal beat classification using pre-recorded data of the MIT-BIH Arrhythmia and 


MIT-BIH Supraventricular Arrhythmia databases. 
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Chapter 3 


Methodology 


In the present chapter, the methodology to find the best combination of features to train the network and 
also the training algorithm is showed. These consists of a theoretical part, in which the ECG signal was 
deeply studied to find parameters characterizing the samples. Then, an experimental part was performed 
in order to select the most suitable training method for the neural network. Once the comparisons led to 
a balanced solution between complexity and overall accuracy. The flow diagram 3.1 shows the step by step 
actions proposed in this work: After data gathering and normalization, the first task is to estimate the 
Bayes error to decide whether to continuing collecting data or extracting the features. Once the features are 
collected, neural network is train by two processed. For the first, non parametric process, particular samples 
are used to design the classifier which consist of the second, parametric process. Finally, the network is 


tested and goes back to gather data. 


3.1 Data Manipulation 


From the proposed methodology, data gathering is the first on the list. In this part of the study there is not 
acquisition of data or instrumentation of patients. In order to manipulate data, a database was obtained 
from the Physionet Bank. Once the data is available, the samples have to be manipulated by the usage of a 


toolbox. 
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Figure 3.1: Proposed Methodology: Flow Diagram 


MIT-BIH Database 


The ECG data used to perform pattern recognition were taken from the MIT-BIH DB: The Massachusetts 
Institute of Technology-Beth Israel Hospital Arrhythmia Database with 48 records, 30 minutes each [13] 
of two channel ambulatory ECG signals. The recordings were digitized at 360 samples per second with 
11-bit resolution over a 10 mV range. In this study, the two channels of ECG signals derived from lead 
I and II were used.Four types of ECG patterns from the database were chosen, which are respectively 
annotated as the normal sinus rhythm (N), premature ventricular contraction (PVC), atrial premature beat 
(A) and left bundle branch block beat (L). Figure3.2 presents the four different HR patterns. After collecting 
the data, cardiologists analyzed the records and then wrote computer-readable reference annotations for 
each beat (110,000 in all the database). Three types of file were released for each record: the reference 
annotator (.atr), the signals (.dat) and the header (.hea). To access and process tthese physiologic signals, 


MATLAB /Octave/Java functions were created and saved as a WaveForm DataBase (WFDB ) Toolbox. 
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Figure 3.2: Different HR Patterns extracted from the ECG cases from MIT-BIH Arrhythmia Database 
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ECG lead configuration 


In most records, the upper signal is a modified limb lead II (MLIT), obtained by placing the electrodes on 
the chest. The lower signal is usually a modified lead V1 (occasionally V2 or V5, and in one instance V4); 
as for the upper signal, the electrodes are also placed on the chest. This configuration is routinely used by 
the BIH Arrhythmia Laboratory. Normal QRS complexes are usually prominent in the upper signal. The 
lead axis for the lower signal may be nearly orthogonal to the mean cardiac electrical axis, however (i.e., 
normal beats are usually biphasic and may be nearly isoelectric). Thus normal beats are frequently difficult 
to discern in the lower signal, although ectopic beats will often be more prominent (see, for example, record 
106). A notable exception is record 114, for which the signals were reversed. Since this happens occasionally 
in clinical practice, arrhythmia detectors should be equipped to deal with this situation. In records 102 and 
104, it was not possible to use modified lead II because of surgical dressings on the patients; modified lead 


V5 was used for the upper signal in these records. 


3.1.1 WFDB Toolbox 


The first beta version (0.9.6.1) was released on February 10th, 2014. The main characteristics of the toolbox 


is that: 


e 31 M-functions 

e Accesses to over 3 TB of physiological signals 

e Interface to LightWave (annotation visualization tool) 

e Multicore processing at Java Virtual Machine (JVM) level (alpha). 


e Unit tested on 64-bit Windows, Mac and Linux. 


The structure of the WFDB data is showed in Figure 3.3. The functions in the toolbox calls a JAVA API. 


After accesing the JVM, it creates and launch a system process to parse, cast, read and write I/O streams. 


The whole process is explained in the Figure 3.4 with the command rdsamp. 
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Figure 3.3: Structure of the WFDB Toolbox 


stdin, stdout, stderr 


Figure 3.4: How a WFDB function (rdsamp) is built. 
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3.1.2 Annotation Analysis 


Annotations are viewed and selected using mnemonic codes rather than the numeric anntyp codes used 
within the annotation files. Annotations can be read by C, C++, and Fortran applications using getann, 
and they can be written using putann, functions defined within the WFDB library. Matlab and Octave 
programs can read and write annotations using these m-files. Alternatively, annotations can be read by 
scripting language applications using rdann, and they can be written using wrann, programs belonging to 
the WFDB Software Package. Details on the format of annotation files are available here. The symbols used 


in plots are describes in Table 3.1. 


Symbol Meaning 


Normal beat 
Left bundle branch block beat 
Right bundle branch block beat 
Atrial premature beat 
Aberrated atrial premature beat 
Nodal (junctional) premature beat 
Supraventricular premature beat 
Premature ventricular contraction 
Fusion of ventricular and normal beat 
Start of ventricular flutter /fibrillation 
Ventricular flutter wave 
End of ventricular flutter/fibrillation 
Atrial escape beat 
Nodal (junctional) escape beat 
Ventricular escape beat 
Paced beat 
Fusion of paced and normal beat 
Non-conducted P-wave (blocked APB) 
Unclassifiable beat 
Isolated QRS-like artifact 
Premature Ventricular Contraction 
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Table 3.1: Annotations used in the Database 
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3.2 ECG Processing 


3.2.1 Pre-processing 
Correction of the signal 
Pre-treatment is to make a clean signal without artifacts. Prior to the removal of noises, it is first necessary 


to correct the signal, this correction is performed by a median filter to correct the base line and reduce the 


signal to a standard benchmark. This operation is shown in Figure 3.5 and by the equations 3.1, 3.2 and 


3.3. 
Yecg2 (N) = Yecgiln) * fi(n) (3.1) 
Yeeg3 (n) = Yecg2 (n) * fo (n) (3.2) 
Yoorrectea (n) = Yecgi (n) a Yecg3 (n) (3.3) 


The median filter is an effective tool to correct the signal without loss of information. This is a filter Classical 
non-linear, its principle is: For each sample of the signal, the class adjacent samples within the window by 
increasing values. We take the median of classified samples and assigns it to the sample current. The 
correction process is developed to remove from the original signal Yecgı (n) the P wave and the QRS complex 
by a first filter median f,(n) of a width of 150 milliseconds, then to remove the T-wave from the Yecg2 signal 
(n) from the first filter through a second median filter a width of 50 milliseconds. Y..g3(n) is the resulting 
signal from the second operation therefore contains as the baseline. The corrected signal Yeorrectea(N) is then 
obtained by a subtraction between the original signal Yecgı (n) and Yeceg3(n) resulting from the two filtering 


operations. 
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Figure 3.5: Effect of the median filter in ECG signal 


Noise Filtering 


The performance of artificial neural network to recognize ECG patterns may be lowered by noise corrupted 
ECG signals. Even though the neural network has some degree of fault tolerance, it is desirable that clean 
ECG signals are provided. Signal is preprocessed , if the sampling frequency is higher then it is downsampled 
and if it is lower upsampled to make the sampling frequency 200 Hz. ECG preprocessing is to eliminate the 
baseline wander and high frequency interference, which can be done by applying a band-pass filtering to the 
ECG signals. Here we used an all integer coefficient digital band-pass filter, where the pass band of the filter 
is from 5-15Hz. Three different type of filters were applied, Gaussian Finite Impulse Response (FIR) and 
Butterworth, shown results can be seen in Figure 3.6 and the corresponding code in appendix A. Matlab 


code of each part of Pre-processing is shown in C. 


3.2.2 Time Series Analysis 


Derivative 


The filtered signal is derivated using a derivating filter to high light the QRS complex (Figure 3.7). 
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Figure 3.6: Filtered signal 115.hea of the PhysioNet MIT-BIH arrhythmia database with QRS markers. 


Squaring 


Signal is squared nonlinearly to enhance the dominant peaks. The result is shown in3.8. 


Moving average 


Signal is averaged with a moving window to get rid of noise (0.150 seconds length). Note : a minimum 
distance of 40 samples is considered between each R wave since in physiological point of view no RR wave 


can occur in less than 200 msec distance. 


R-peak Location 


R peak detection in QRS complexes is well searched since it is considered to be the basis of ECG wave 
segmentation, and several algorithms have been developed to locate the R point, some of which acquired a 
low error rate less than 0.5 % [3]. Depending on the sampling frequency of the signal, the filtering options 
can changed to best match the characteristics of the ecg signal and the R location, used in the present work, 


to find the QRS complexes. The code used to find peaks is described in appendix B. The results of this 
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Figure 3.7: Derivative filtered signal 115.hea of the PhysioNet MIT-BIH arrhythmia database 
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Figure 3.8: ECG vs Squared ECG 
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findings are shown in Figure 3.9. Initially FFT is applied on the ECG signal using equation 3.4. Then low 


frequency components were removed. On the resultant signal inverse FFT is applied given by 3.5. Now the 


signal is filtered for detecting the R peaks. 
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Figure 3.9: R-peak Location 


Peak Detection 


The R-point, as calculated above, is found to have reasonable accuracy. Q wave is defined as A negative 
wave at the onset of the QRS complex and the valley [minimum] is defined as Q point. Thus to position 
the Q point, it is positioned as the local minimum in a short (approx. 0.05 sec) window around the left of 


R-point estimated in 3.6. 


Uogys and W ogy asat (3.6) 
dg dx 


S Peak Detection 


The S-point is first approximated as the point where the slope has its first negative to positive zero crossing 


after the R-point found in . It was found that it could be placed more accurately as the local minimum in a 
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0.05 sec window before the above approximation (See . 


d 
TY co y x<a(s) and 
dx 


dy 


£ 


>0 Y 2>x(s) (3.7) 


This is to be true for the s window length of 0.05 sec. 


3.2.3 Window Finding 


To find the QRS complex window, one have to extract the most important features of the signal in order 
to be able to reconstruct it. This section describes the QRS complex detection, rhythm analysis and the 


resulting function for the window finding. 


QRS Complex Detection 


The first step in almost any ECG analysis procedure is the detection of QRS complexes. QRS complexes 
are the most distinct parts of the ECG and compared to other signal parts they can be detected rather 


easily. All the time series analysis has the objective of finding specific features to find QRS complexes: 


1. High amplitude signals 


2. High derivatives 


A simple threshold based QRS detector is described below. Examples of more sophisticated algorithms are 
described in [26] and [31]. The more sensitive and specific parameters of these two is the derivate. Therefore, 
calculate the first derivative of the signal. Since there are ECGs, where the R peak is more distinct, and 
others where the S peak dominates, and since sometimes the ascending and sometimes the descending part 
of these peaks is more narrow, both positive and negative extrema of the first derivative of the signal can 
be used for detecting QRS complexes. Therefore, the absolute value of the first derivative of the signal is 


often used for QRS detection. The resulting signal can be summarized by the equation 3.8. 


Sigaet = abs(dif f(filt(sigorg))) (3.8) 


To the resulting signal a threshold can be applied, assuming that all over threshold values are situated 


within a QRS complex. This threshold method is a very simple kind of QRS detector. Alternatively, the 
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Figure 3.10: Absolute value of the first derivative of the filtered signal 115.hea (name of data-set) of the 
PhysioNet MIT-BIH arrhythmia database 


Pan-Tompkins [26] algorithm could be used: Start at the first sample. Compare the sample to a threshold. 
If it is lower than the threshold, take the next sample. If the samples amplitude is higher than the threshold, 
interpret it as a QRS complex. For the samples right after such a super-threshold sample use an increased 
threshold that decreases exponentially to the original threshold (refractory period). Take the sample right 


after the super-threshold sample and compare it to the (increased) threshold. Usually, more than one sample 
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Figure 3.11: Absolute value of the first derivative of the filtered signal 115.hea of the PhysioNet MIT-BIH 
arrhythmia database. Values ¿ 0.1 are marked with red asterix. 


within each QRS complex is super-threshold as shown in Figure 3.11. In order to get one single sample per 
QRS complex, the original filtered (but not derivated) signal can be used. Super-threshold samples that 
are closer than e.g. 20 ms to one another can be grouped in super-threshold regions. Within each of these 


regions the maximum of the original signal can be interpreted as the R peak of the signal. 
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Rhythm Analysis 


Soon as the QRS complexes are available, the RR intervals of the ECG can be calculated simply by calculating 


difference in between the time points where two subsequent R peaks were detected (Equation 3.9). 


RRi = tr... — tr, or RR = dif f (tr) om 


Figure 3.12 shows an example of RR intervals plotted over beat numbers and is also called Tachogram. Note: 
this figure has been generated without manual corrections of the results from the QRS detection algorithm 
and contains therefore sudden changes because of wrongly assigned QRS complexes. Although these RR 
intervals are varying from beat to beat (Heart Rate Variability, HRV), the variations are normally quite small 
(only a few percent). Therefore, ventricular rhythm abnormalities can be found by finding RR intervals that 


vary more than just a few percent from the preceding RR intervals, e.g. more than 15 percent. 


abnormalbeatindices = find(abs(dif f(RR)) > RR/15) (3.10) 


Unfortunately, QRS detection seldom works with 100% sensitivity and 100% specifity usually some QRS 
complexes are missed and some artefacts are interpreted as QRS complexes. Therefore, manual inspection of 
the detected QRS complexes and manual corrections may be indicated before detecting rhythm abnormalities 


according to the equation above. 
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Figure 3.12: Tachogram, i.e. RR-Intervals as a function of time of signal 115.hea of the PhysioNet MIT-BIH 
arrhythmia database. 
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Window finding algorithm 


A function for window finding was implemented following the steps described in figure 3.13. The inputs and 
outputs of the proposed system are shown and also the complete algorithm can be found in appendix C. 


Also, the resulting window can be seen in figure 3.14 
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Figure 3.13: Window Finding proposed system: Inputs and Outputs 
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Figure 3.14: Found window for 110.hea from MIT-BIH database 


Heart Rate Calculation 


Feature as RR interval (used medically as an indicator of Ventricular Heart Rate) metrics are generated 
from the deflection positional information. To determine RR interval two R, peaks in consecutive beats is 


calculated and their difference is computed. Heart Rate is 60/RR interval beats per minute (See 3.11). 


60 : 
Rate = Pini beatsperminute (3.11) 


3.2.4 Compresion of the ECG signal by linear prediction 


The basic idea of feature extraction method characteristics of the ECG signal is the data can be obtained by 


a linear combination past data of this signal. The sample predict is calculated from samples passed through 
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the relationship shown in equation 3.12. 
i p 
Yn, =X ai Yn- (3.12) 
i=1 


Where a; is the ith coefficient of the prediction Linear (LPC: Linear Prediction Coding) and p its order. 
To reduce the size of the input vector neural network, and using the work of Lin [19] which showed that 
approximately equal prediction two (P = 2) was sufficient for signal analysis ECG, we have implemented 
this compression method. The QRS complex is packed in two coefficients (a, and a2) and the P wave (az 
and a4). A general survey of work on ECG signal classification [1],[13] and [30] led us to choose a better 
representation of the complex QRS and P wave, namely: The coefficients (a; and a2) of the resulting applying 
the LPC plus RR interval and QRS width are the elements of vector representing the QRS complex of each 
cardiac cycle. The information in each cycle are represented by a vector of four elements the QRS complex, 
five components for the P wave The coefficients are the LPC morphology of the two waves, while the other 


elements are connected temporal characteristics the designated cycle. 


3.3 Abnormalities 


Abnormalities originating from the ventricles do usually have abnormal QRS morphologies, since the conduc- 
tion spreads over the ventricles in a way different to normal beats. Therefore, by comparing the morphology 
of a beat with abnormal RR interval to that of a normal beat, one can decide, whether the rhythm abnor- 
mality’s origin is within the ventricles (ventricular rhythm abnormality) or in a tissue normally excited prior 
to the ventricles (supra-ventricular rhythm abnormality). 

When comparing the morphology of a beat to that of a normal beat, the normal beat’s morphology needs to 
be known. One possibility for calculating a normal beat’s morphology is to determine the mean morphology 
of all heart beats. As long as most beats within the ECG are normal and only little beats are abnormal, the 
average will correspond to the normal beat’s morphology. Better results will be achieved when averaging 
only beats with normal RR intervals but still the result will only be correct as long as the excitation of most 
beats with normal RR interval spreads over the ventricles in a normal way. One possibility for checking the 
similarity of two heart beats is the cross-correlation. The cross-correlation in between two signals can be 


calculated with the following equation 3.13. 
N 
cc(n) = 5 beatı (n) x beata (n) (3.13) 
i=—N 
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where N is the number of samples per beat. 

For each n its value varies in between -1 and 1. Zero corresponds to no correlation in between the signals. 
-1 results from signals that are mirrored by the x-axis. 1 means that the signals are identical. cc(n) 
corresponds to the correlation in between the beats when displacing beatz n samples to the right. Therefore, 
cc(0) describes the similarity of the two beats without adjusting their position along the time axis and can 
be used for deciding, how similar the morphology of two heart beats is. The correlation factor in between 
two similar beats is usually higher than 0.95. 

Figure 3.15 shows a 2 channel ECG recording with a ventricular extrasystole. The RR interval is reduced 
and the morphology is different to that of the QRS complexes with regular RR interval. No P wave can be 


seen. 
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Figure 3.15: Signal 121 of the arrhythmia database with a ventricular extrasystole in the middle of the figure. 


3.3.1 Arrhythmia Detection 


The abnormalities in ECG morphology lead to find a particular case of arrhythmia, in this study three types 
of arrhythmia are being studied for an automatic detection: PVC, A, L (See table 3.1). This is to be done by 
the usage of a classifier, which receives the features as an input and finally decides (by classification) which 
type of arrhythmia the ECG is coming from. A diagram describing the proposed system is shown in figure 
3.16. From section 2.2, we proposed three types of arrhythmia described in table 3.2. With these features 


various Cardiovascular Arrhythmia are detected as Right bundle branch block. Right bundle branch block is 
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Arrhythmia ECG signal 


PVC 


N 


Table 3.2: Arrhythmia studied in this work 


a delay or block of conduction within the right bundle branch. A delay of conductionmanifests as incomplete 
right bundle branch block. A QRS duration greater than 0.14 sec results in right bundle branch block. 
Bradycardia occurs on resting heart rate of under 60 beats per minute, though it is seldom symptomatic 
until the rate drops below 50 beat/min. Tachycardia refers to rapid beating of the heart as a heart rate 
greater than 100 beats per minute in adult. Diseases were predicted from these features derived as according 
to medical science for Tachycardia heart rate > 100 bpm(beats per minute). Normal Sinus Rhythm is 60 to 
100 bpm. Ventricular tachycardia for heart rate is from 101 to 250 bpm and QRS width > 0.12 sec. Normal 
QRS width is 0.04 - 0.10 sec. Incomplete Bundle Branch Block for QRS width between 0.10 sec and 0.12 


Raw ECG 
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Figure 3.16: Classification system proposed. 
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sec and Bundle Branch Block for QRS width > 0.12 sec. 


The questions to be answered for this study are: How is the classifier going to decide? and What are 


going to be the features?. Sections 3.5 and 3.4 respectively. 


3.4 Feature Selection 


The first step was to analyze the signal in time by the extraction of the most important characteristics of the 
signal to generate an adequately accurate QRS complexes window. Finding features lead to represent the 
original signals with lower dimension, which turns to form the input vector space. The ECG patterns 
to be recognized exhibit dissimilarity both in heartbeat timing and ECG morphology. Study done by [13] is 
taken as an example for feature selection, in this particular case, the number of features used was a total of 
29, consisting of a ECG sequence pf 50ms before and 100ms after the R point, QRS morphology with no P 
and T wave, two R-R intervals and R peak location. In the present research different tests were done in order 
to find the best of the feature combinations. To propose the features, was a work of study of the signals an 
intuition based on the knowledge of the behavior of the abnormalities. Table 3.3 shows the observed features 
over different examples used for this study and plotted from the mitdb database. For each of the arrhythmia 
and also for normal beats, the following characteristics were calculated and checked how they vary in order 


to realize how important they are to lead a classifier to a right decision: 


e P-peak location (X axis) 

e R-peak location (X axis) 

e Previous RR interval (X axis) 
e Next RR interval (X axis) 

e Q-S duration (X axis) 

e Q-T duration (X axis) 

e P-R duration (X axis) 


e Heart Rate Variability (X axis) 


42 


Arrhythmia | ECG signal 
A 
e Rate: Varying from stop to rapid (350 to 650 bpm) 
e Rhythm: Irregular 
e P wave: Fibrillatory (fine to course) 
e PR Interval: None 
e QRS complex < 0.12 seconds; wide and bizarre and followed 
by a compensatory pause 
PVC 
e Rate: Depends on underlying rhythm 
e Rhythm: Regular or Irregular 
e P wave: None 
e PR Interval: None 
e QRS complex 0.12 seconds; wide and bizarre and followed 
by a compensatory pause 
0.15 seconds | | | 
No R 
L 
e Rate: Depends on underlying rhythm 
e Rhythm: Regular 
e R wave: None 
e S wave: Deep 
e QRS complex: 0.12 seconds 


Table 3.3: Characteristics of the arrhythmia seen from the V1 Lead Instrumentation 


e P-peak (Y-axis) 


43 


Q-peak (Y-axis) 


S-peak (Y-axis) 


T-peak (Y-axis) 


U-peak (Y-axis) 


ST segment (2) 


In the present study, and at this stage 14 features have been proposed. However, the purpose of this 
study is to reduce even more the number of features, so that different offline and online experiments will 
be achieved in order to test the best of the combinations. It is also important to understand not only the 


characteristics for arrhythmia ECG signals but also for normal beats. 


3.4.1 Characteristics of a Normal ECG beat 


Figure 3.17 shows an expected ECG signal from a normal heart beat. Similarly, table 3.4 shows the featured 
values for this waves. 
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Figure 3.17: Normal Beat ECG 


3.4.2 Arrhythmia Prediction Using ANN 


The artificial neural networks play a significant role in the field of artificial intelligence. In spite of the advent 


of computers, a system which can imitate the human brain is of a great demand. Artificial Neural Network 
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Parameter Normal Value Remarks 


Heart Rate | 60-100 bpm Atrial (A) and Ventricular fibrillation (V) causes stroke and heart 
quivers leading to sudden death of person. 
P wave Amplitude: 0.25 + 0.05mV and | Electrical activity associated with the contraction of atria. 


Interval 110 + 20ms 
QRS complex | Amplitude: 1.6040.5mV and In- | Associated with ventricular contraction. 
terval 100 + 20ms 
R wave Amplitude: 1.6 u 0.5 mV — 


Q wave Amplitude: 0.25 times the R | — 
wave. 
T wave Amplitude: -0.5 mv and Interval: | — 
160 ms 
PQ interval | Interval: 120 - 200 ms Time taken by the electrical signal to travel from atria to ventricle. 
RR interval | Interval: 120 - 200 ms Time taken by the electrical signal to travel from atria to ventricle. 


Table 3.4: Parameters for Normal Heart Beat ECG 


is the technique where we try to copy the working of human brain. It has a very significant role in the 
field of artificial intelligence. An ANN comprises of interconnection of artificial neurons which follows the 
function of biological neurons and are basic building blocks of the network. They learn from the data fed 
to them and keep on decreasing the error during training time and once trained properly, their results are 
very much same to the results required from them [16] thus referred to as universal approximators [11]. The 
error is decreased by modifying weights of individual neurons. The most popular neural networks used by 
researchers are the multilayer feed forward neural network trained by the back propagation algorithm [16], 
[15]. There are different kinds of neural networks classified according to operations they perform or the way 
of interconnection of neurons. The ANNs are capable of learning the desired mapping between the inputs 
and outputs signals of the system without knowing the exact mathematical model of the system. Since the 


ANNs do not use the mathematical model of the system, they are excellent estimators in non linear systems. 


3.5 Classifier 


To perform a classification is to set identification of relationships between objects and between these objects 
and their parameters. From nearby or dissimilarities, it is to build a partition of the set of objects in a set of 
classes as homogeneous as possible. The work of research for the classification of cardiac signal date back to 
1957 when the first approaches were mainly based on the use of the method decision trees and the statistical 
methods. Thereafter, new approaches have been developed, such as the syntactic approach, expert systems, 
fuzzy logic and network neurons. To this end, we propose to design a classifier of cardiac arrhythmias, classes 


output are: 
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e Normal (N) 


Premature Ventricular Contraction (PVC) 


Atrial premature beat (A) 


Left bundle branch block beat (L) 


To solve a classification problem, we have a range of methods to be used by knowing their capabilities and 
especially their limitations. Throughout the experiments (theoretical and practical) carried out by several 
researchers [21], neural networks lead to good results. Properly sized, they get always the best classification 
rates. A big number of parameters such as variations of heart rate, width and morphology QRS and P 
and T waves of the signal are necessary, and are used to assess whether the patient presence or absence of 


anomalies. 


3.5.1 Software Implementation of the classifier 


The purpose of the software implementation of the classifier is to determine the size and parameters of the 


neural network, namely: 


e Number of layers and number of neurons for each layer. 
e Error reached. 


e Number of iteration. 


These parameters are used to have the best network performance. In this section, a discussion on the 


conditions and method of learning, programming and dimensioning of the network is performed. 


Learning 


Learning seeks to discover relationships and hidden patterns in data collections. A learning system takes as 
input a set of examples (cardiac cycles) seeking for a classification. In our case, the system trying to learn 
the forms of arrhythmia heart from the examples provided at its entrance. These examples are presented as 


parameters (elements of the input vector) which and represent the temporal characteristics morphological 
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cardiac cycles classified according to arrhythmia which they correspond. The Back-Propagation learning 


algorithm in this case ensures that the generated classes allow discrimination to better examples entry. 


Back Propagation Algorithm 


Supervised neural networks require Learning by trial/error, starting examples given in the form of pairs 
(problem, solution to this problem). The numerical method used for learning and adopted in this study is 
called gradient descent backpropagation. The idea of this optimization algorithm is that the gradient of a 
function always points in the direction where this function is growing the most. Vice versa, If we multiply 
by the scalar -1 gradient, we get the direction of the greatest decrease function. The major difficulty is 
choosing the no learning. Indeed, if we are going to choose too small, we will converge to a minimum Local 
very slowly and if we choose too big, we do not converge, because we are going jumping from one side to 
another. Several methods and optimization algorithms have been proposed to solve this [12-14] problem. 
Each method has advantages and drawbacks. Of course, no method led fail to global minimum. it should 
therefore be placed under conditions in which the minimum premises are as few as possible. We chose the 
Levenberg-Marquardt method that is to change the settings according to the relationship showed in equation 


3.14. 


Ok = 0x1 — [hk-1 + Arp-11] (3.14) 


Where I: identity matrix. This method is particularly clever because it adapts itself to the form of the cost 


function. It carries a compromise between management gradient and direction by the method of Newton. 


3.5.2 Training Data 


Signals from 48 records MIT-BIH Database (See 3.1) are used for the development and evaluation of classifier 
performance. Each registration of this database is accompanied by an annotation file in which each ECG 
beat (cardiac cycle) was identified by experts (cardiologists). These labels (annotations, see table 3.1) are 
used to learning classifier and evaluation of these performance during the test phase. Since this study aims to 
identify the premature ventricular (PVC), atrial fibrillation (A), left bundle branch block (L) some database 
records of PVC-free, A-free, L-free data were excluded from the study, leaving 35 interest records. Heart 


cycles used for learning different pathological cases were randomly selected from several recordings. The 
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advantage is to build on a database including different forms of each cardiac arrhythmia. 


3.5.3 Network Sizing 


There is no theoretical results, or even satisfactory empirical rules, which allow properly size a neural network 
function problem to solve. Designing a multilayer network is experimentally, the difficulty usually arises when 
selecting the number of intermediate layers and number of neurons in each. For the input layer, it contains 
as many neurons as the number of elements of the input vector. Similarly, for the output layer, it contains 
as many neurons as the number of classes to be discriminated. Therefore, the input layer consists of nine 
neurons and output layer consists of six neurons. To determine the number of neurons hidden layer, we 


proceeded in the way as follows: 


1. Prepare the cardiac cycles corresponding to arrhythmia chosen for learning. 


2. Create the neural network whose number of the input layer (of respectively output) is set according to 


the input vector (respectively in output level). 
3. Appoint an arbitrary number of neurons in the hidden layer. 
4. Attach a mistake and very low value the number of iterations. 


5. Start learning. 


As the process diverges, number of neurons in the hidden layer increases, if the next one is too high, the 
number of intermediary layers increases. Unfortunately, there is no super neural architecture, that is the 
reason why it would be sufficient to adapt the weight problem in order to get the best performance possible. 
In general, the size and type of connections of a network influence its performance. Note however that this 
architecture found is not unique, and there are often several architectures to obtain the same performance. 
In our case, the optimal architecture will be the one with the smaller, or more specifically one that will 


require fewer calculations to perform the requested function. 


3.5.4 Network Architecture 


To solve the problem of arrhythmia identifying, a Feed-Forward Network with one input layer, one hidden 


layer and one output layer is proposed as shown in the Figure 3.18. The input layer consist of five neurons 
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Figure 3.18: Network Architecture 


with the transfer function of tan-sigmoid, the hidden layer consist of four neurons with the transfer function 
of log-sigmoid whereas the output layer consist of four neurons with the linear transfer function. The inputs 
given to the input layer are R-peak location in x and y axis, R-R previous and next interval, QRS duration, 
Heart Rate, U peak and ST slope whereas the outputs obtained are the presence of arrhythmia namely 
Atrial Premature Beat (A), Premature Ventricular Contraction (PVC) and Left Bundle Branch Block Beat 


(L). Where, category 1 is A, category 2 is PVC and category 3 is L as shown in table 3.5. Hence the output 


Beat Symbol Type of Beat Category 
A Atrial Premature Beat 1 
PVC Premature Ventricular Contraction 2 
L Left Bundle Branch Block Beat 3 
N Normal Beat. 4 


Table 3.5: Categories by Beat type 


from the network will be governed by the following equations 3.15, 3.16 and 3.17. 


4 
Tia = tansig(O Wiaj X P;) + bia) (3.15) 
j=1 
5 
Toa = logsig( O ` Waj X T1;) + boa) (3.16) 
j=l 
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4 
Tsa = Da = purelin((S~ W3aj X Toj) + bsa) (3.17) 


j=1 


h th 


Where Tmn is output from nt” neuron of m*” layer, Wpor is weight for j*” input for q*” neuron of pt” layer and 
pq 8 J 


h neuron of mt” layer. The network is trained by Levenberg Marquardt Back-Propagation 


bmn is bias for nt 
algorithm which updates the weights and biases according to Levenberg Marquardt up to the mean squared 
error of 0.001 which was attained in one hundred and seventy eight epochs with 20 training dataset. The 
training graph is shown in Figure 3.19. As soon as the training procedure is over, the neural network gives 
almost the same output pattern for the same or nearby values of input as visible by the mean squared error 


in Figure 3.19. This tendency of the neural networks which approximates the output for new input data is 


the reason for which they are used as intelligent systems. 


Best Training Performance is 0.000927112 at Epoch 200, Goal is 0.001 


Mean Squared Error (mse) 


200 Epochs 


Figure 3.19: Training of the Network 


Training Algorithm Design 


Our choice was focused on the backpropagation algorithm, a program was written in Matlab for the imple- 
mentation (Apendix D). The activation function is selected by the sigmoid function (Figure 3.20). A sigmoid 
function is defined by f(x) given by equation 3.18. 


e*z _ 1 


= 1 
Oka +1 (3 8) 


f(x) 
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Figure 3.20: Sigmoid Activation Function and its firts derivative 


Indeed, it goes to a as x approaches +00 and goes to —a when x approaches —co. The parameter a then 
regulates the saturation value; the parameter k is used to control the slope of the curve any point outside 
saturation. In the presented model, functions logsig and tansig are implemented as Matlab activation function 
(transition) in the network neurons. These functions are continuous and sigmoid differentiable, they are also 
not decreasing. We call quasilinear functions because they are almost in linear the range of interest. As 
shown in the figure 3.20, sigmoid function f(x) and its first derivative f'(x) are both continuous functions. 


This is useful when handling mathematical functions. 


3.6 Arrhythmia Recognition Module for Mobile Devices 


Arrhythmia is something that may happen to people at any time or place or while performing any type 
of activities, leading to fatal outcome in number of cases. Tracking vital functions and getting a response 
from a system can help people stabilize their organism, make some precautions when unexpected behavior 
arises or take some necessary actions while waiting for medical help. Most of the proposed systems [23] are 
not interested in the type of the arrhythmia - they only want to raise an alert signal when something is 
wrong with the heart rhythm and behavior. The proposed solution in this study makes diagnosis for the 
users condition, aimed to advice someone to give first aid assistance while waiting for medical help. It is 


not meant to trigger a whole system, like calling an ambulance, upon raising an arrhythmia signal 3.21. It 
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Figure 3.21: An overall look of the system 


is a great concern for this system not only to work with satisfactory reliability but also with low resource 
consumption, so that other applications, that are commonly used while performing this kind of activity, 
could work in parallel: Internet radio, MP3 players, GPS navigation or other applications, or a combination 


of them. 


Proposed Overall System 


As shown in figure 3.21, the final intended system is comprised of ECG transmitter, vital signals sensors and 


Mobile ECG Application. The mobile application has two main modules: 


e ECG Processing Module 


e Detector Module 
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The system consists of Feature Extractor and Arrhythmia Recognition implemented as a mobile application 
and is intended to be installed on smartphones. The intenyion is to raise alert signal and display a short 
text message, after recognizing a possible critical state of the patient. The signal comes from the ECG 
transmitter and is then passed trough mobile phone, where pulses are extracted and tested against trained 
arrhythmia recognition machine. If an arrhythmia is recognized, an alert signal is raised and an emergency 
call is launched. Along with the data from the vital signals sensors and trough interaction with modules the 


user gets valuable information of how to engage the heartbeat disorder. 
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Chapter 4 


Results 


During the development of the present work, offline, online and mobile platform implementation were per- 


formed, obtaining results and conclusions described as follows. 


4.1 Window Extractor 


A function for window extraction was implemented in order to send only the important characteristics of 
the signal to the mobile platform with the neural network algorithm implemented on it. The results of 


the window extractor can be seen in figure 4.1 and the algorithm is to be consulted on C. Windows, being 


e S Buffer 
| 
) 
INPUT 
i i Signal Cleaning 
| l! | & Formatting 
OUTPUT 


Pulse Extractor 


Figure 4.1: Stream preparation for window extraction 


either regular or irregular (arrhythmia), consist of number of samples ranging from 200 to 500 samples per 
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window. We included 35 different patients (recordings) in the process of obtaining pulses. We used 25 of 
them for obtaining arrhythmia pulses, and 10 of them for obtaining regular pulses. The pulse extraction was 
performed automatically with a light version of the Pan & Tompkins algorithm (without ECG cleaning and 
formatting) and some manual calibrations of the parameters of the algorithm. This algorithm basically finds 
a piece of the stream with defined length, where each of the slopes between points at a certain distance is in 
predefined bounds. This way we detect the separation of the pulses. After-wards we detect the closest slope 
above some value on the left and on the right of the separation line. The left point marks the end of the 


previous pulse, and the right one the beginning of the next pulse as shown in Figure 4.2. When extracting 


Figure 4.2: Finding the separation line, the end point and the start point 


regular windows we must make scaling of the both X and Y axis. Scaling of the Y axis eliminates the vertical 
offset of the pulse (people do not have same potential when heart is in relaxed mode). This will eliminate 
the false features that may occur because of that offset: we are interested in the pattern of the regular 
pulse, not the offset. Arrhythmia can be scaled only by the X axis. Since we do not make classification 
of the arrhythmia by their type, the difference in the offset is not important. The window contains all the 
important information about the signal. The maximum number of components of this characteristic vector 
was proposed to content 15 elements (See section 3.4). To reduce the length of this vector is a target of 
this work. Window extraction for the worst case ECG signal which involves lots of noise with it and high 
variations is taken to prove the efficiency of the system. The various steps involved in Section 3.2 are shown 
in Figure 4.3 for the worst case which clearly shows that the algorithm is very good in handling even the 


worst cases especially because of its two pass filter. 
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Figure 4.3: Window Extraction Result 


4.2 Implementation of the Neural Network 


Since the process is divided on two levels, first identifying the features and second predicting the disease 
from the features identified. Hence, the accuracy of the results depends on accuracy of both levels. So the 


results on both levels have been depicted below. 


4.2.1 Feature Reduction 


Offline implementation lead the best feature combination after a number of experiments which were per- 
formed with the usage of MIT-BIH database. From works such as [20], [7] or [28] the number of features was 
not so important to be reduces since they were trying to compare different types of neural networks with no 
care on complexity increasing. In this study, both, the overall accuracy and the simplicity, were taken into 
account. 15 features were proposed leading to a total of 77,558,760 to be performed. Since we are not able 
to perform such a huge amount of experiments, we performed a feature extraction for the worst case with 
three different neural network models.4.1. The selected features were, R location on the x-axis and y-axis, 
RR previous and next interval, QS interval, HRV, U location and the ST transition slope. Thus, 8 features 


were selected at the end, reducing complexity of the network compared to other studies such us [20]. 
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Feature Trial 1A | Trial 2A | Trial 3A | Trial 1B | Trial 2B | Trial 3B 

P X X - - - - 

R X - X X X X 

Previous RR interval X X X X X X 

Next RR interval X X X X X - 

Q-S duration X - - X - X 

Q-T duration - - - - - - 

P-R duration - - - - - - 

HRV X X - X - - 

R-peak X X - X - X 

Q-peak - - - - - - 

S-peak - - - - - - 

T-peak - - - - - - 

U-peak X X X X X - 

ST segment (5%) X E X x X X 
Overall Accuracy SOM 90.14% 76.21% 83.12% 91.12% 89.04% 90.31% 
Overall Accuracy LVQ 88.12% 74.54% 87.10% 90.03% 88.22% 91.32% 
Overall Accuracy BP 89.90% 75.02% 80.76% 91.78% 87.64% 90.89% 


Table 4.1: Characteristics of the arrhythmia seen from the V1 Lead Instrumentation 


4.2.2 Neural Network Complexity 


Similarly, several combinations were tried with different number of neurons in the hidden layer. The results 


of these changes are depicted in the Table 4.2. The decremented performance on increasing the number of 


Number of Neurons in | Accuracy 
Hidden Layer 

2 89.23% 
3 90.26% 
4 91.70% 
5 91.05% 
6 90.67% 
7 90.67% 


Table 4.2: Effect of hidden layer neurons on accuracy. 


neurons could be due to over-training. On the basis of above readings, the network with five neurons in the 
input layer, four neurons in the hidden layer and six neurons in the output layer was selected. The network 
so formed was then tested on 35 dataset with an accuracy of 98.48%. Figure 4.4 shows the predictions and 
the actual disease identification for each of the thirty five dataset. X-axis represents the dataset id whereas 
the Y-axis represents the disease id. A circle represents the actual value (target value) and the star (*) 
represent the value predicted b the network. The results shown above clearly indicate that the system can 


assist a doctor to a very good level about disease indications from the ECG and hence the necessary steps 
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could be taken before the disease strikes the patient. The proposed solution has an accuracy of 91.70%. 
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Figure 4.4: Predicted output against desired output. 


Table 4.3 depicts a comparative study of various techniques used previously and clearly indicates that the 


proposed method competes with overall accuracy of previous ones but it is much simpler than them, due to 


feature reduction. 


Method Accuracy | Elements on Feature Vector Neurons in the Hidden Layer Total Neurons 
SOM [20] 95.50% 29 64 neurons (Competitive Layer) 129 
LVQ [20] 91.50% 29 10 43 
BP [20] 92.50% 29 10 41 
BSS-Fourier [27] | 85.04% 3 6 24 
MOE [17] 94% 15 15 36 
Fhyb-HOSA [6] 96.06% 9 22 36 
DWT-NN [4] 96.79% 18 12 24 
FTNN [22] 92.0% 5 82(Competitive Layer) 214 
ICANN [34] 93.37% 7 6 17 
Proposed 91.70% 8 4 16 


Table 4.3: Comparative results of different methods. 


4.3 Testing of the proposed network 


4.3.1 Automotive/Aerospace Neurophysiology Lab 


Previous offline implementation was tested with same data which one the network was trained. In order to 


test the proposed solution in an online environment, an automotive/aerospace neurophysiology Lab, proposed 


by gTec, was used. The presented laboratory served to collect signals from patients on-line in such a way 


that the proposed algorithm could be tested on real time. 
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Automotive/Aerospace Neurophysiology Lab 
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Figure 4.5: Automotive/Aerospace Neurophysiology Lab 


4.3.2 g.BSanalyze software 


g.BSanalyze is a patented software, free for g Tec hardware buyers. It aims to manipulate all the components 
in a common instrumentation such us the one presented in Section 4.3.1. One can adjust the sampling 
frequencies, filter and do any mathematical operation with the signals since the program is Matlab based 
and contains a number of functions plus the ones given by the already implemented toolboxes. To access to 


the modules, special permission was given. 


Filtering Data 


The data-set 100.dat can be found under ECGAcquisition/Test100. The captured data segments are also 
exported to the MATLAB workspace export115.mat. To filter the input data, the MATLAB fft and ifft 


commands are used (batcha1.m): 
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load ('export100.mat' ); 
seql=data(: ,1); 
seq2=data(: ,2); 
figure(1); clf; hold on 
subplot(2,1,1); 
plot (seq1(1:4000) ,'k') 
subplot (2,1,2); 
plot (seq2(1:4000),'k'); 


nsamp=length (seql); 
NFFT = 2" nextpow2(nsamp); 

f = fs/2*linspace(0,1,NFFT/2+1); 

highpass = 1; % highpass frequency 

lowpass = 10; % lowpass frequency 

% process channel 1 

freq_seql = fft(seql, NFFT) / nsamp; 
wrong.inds = find(f < highpass | f > lowpass); 
wrong.inds = wrong.inds(wrong.inds <= nsamp/2); 
wrong_inds = [wrong.inds, nsamp - wrong_inds]; 
freq_seql (wrong_inds) = 0; 

filt_seql = real(ifft(freq_seql)) * nsamp; 
filt_seql = filt_seq1 (1:nsamp); 

% process channel 2 

freq_seq2 = fft(seq2, NFFT) / nsamp; 

freq-seq2 (wrong-inds) = 0; 

filt_-seq2 = real(ifft (freq_seq2)) * nsamp; 
filt_seq2 = filt_seq2 (1:nsamp); 

figure (1); 

subplot (2,1,1); hold on 

plot (filt_seq1 (1:4000),'r') 

subplot (2,1,2); hold on 

plot (filt_seq2 (1:4000),'r') 
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QRS Complex Detection 


Window Extraction algorithm has been used to identify features from an ECG signal. Using the known 
characteristics of the network, prediction of arrhythmia was achieved. It is found that the system is very 
robust and can identify and predict features even from highly abnormal ECG. This is a big benefit since the 
ECG pattern varies in many factors from person to person. These factors may be height of peaks, width of 
QRS complex, presence or absence of peaks, heart rate etc. QRS complex is the most important section of 
any ECG and once it has been detected, one can use it in other study or system also. High accuracy of the 


system makes it highly reliable and efficient. 


e Start the Complex Detector Beat Analysis to find QRS complexes (See Figure 4.6) 


e Evaluate the QRS complexes Detector by Correlation Analysis 


importa = x 
Lookin: | |, Lesson4 +] «#6 Ck E 


Name | +] Date modified | +| T, =| Size +] Ti X 
1100,hea 1 103,hea 
11.hea 115.hea 
L121,hea 1 J122.hea 


Desktop 
"mm 


Guger 


a 


File name: [115hea +] 
Files of type: [an DAT-Files (*.hea] +] Cancel | 
A 


Figure 4.6: Importing Data in g.BSAnalayzer 


After inserting the proposed code for window extraction, a correlation is run in order to evaluate the algorithm 
presented. The correlation performance analysis from the ECG menu can be seen in Figure 4.7. The Single 
Beat Editor opens with the correlation result. Click on Corr. Class 1 to view normal beats (Figure 4.8): 
Click on all different classes to investigate the different complexes. Class 3 is visualized in Figure 4.9: Use 
the Jumper function in g.BSanalyze to visualize the different QRS complexes corresponding to the different 
classes. Select e.g. QRS 3 to see the complex at second 1518 (Figure 4.10). The correlation classes show the 
averaged beats for each specific QRS number. The Automatic classification assigns new markers with the 


specific pathologies to the QRS complexes. Seen in Figure 4.11. 
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Bikes 


Classify all selected markers according to the signal morphology within a certain time 
window. Markers with similar signal morphology (with a correlation coefficient that is 
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Figure 4.7: Correlation Analysis of the Window Extractor 


Filename: CECGLecture\Lesson4it 


4.3.3 Acquiring ECG Signals 


Car Driving Simulation 


The instrumentation of a patient was performed only to acquire ECG signals for online testing of the network. 
The relaxed ECG signal can be seen in figure 4.12 and when the car driving simulation starts, the ECG can 


be seen 4.13. 


4.3.4 Detection of arrhythmias and abnormalities in the conduction system 


using gTec laboratory 


The first step was to import the data in g.BSanalyze. Scroll through the data-set to find points of interest. 
Write down the time point and a description of the signal/artifact/pathology. Points of interest (time point 
and description). All the experimental data and the diagnose is shown in figure 4.14. The ECG records 
from 4 different types of patterns were obtained from 13 patients (patients numbered as 100, 101, 107, 
108, 109, 110, 111, 112, 113, 114, 115, 116, 117). 176 ECG segments were chosen for each of the pattern, 
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Figure 4.8: Single Beat Editor 


Figure 4.9: Class 3 Correlation 


which produces a dataset of a total number of 704 ECG records, each containing a QRS complex. The 
four patterns are respectively designated as N, A, L, PVC. We used the back propagation neural network, 
that is, the whole dataset was divided into four groups of equal size with equal number of 4patterns in each 
group, and then the proposed neural network was tested by a different data group, while the other three 
is used for training. Before being fed into the neural network models, the ECG signals were preprocessed. 
The ECG preprocessing amd yhe neural network were simulated in Matlab Mathworks with the usage of 
g.BSanalyzer. The performance of the online neural network is evaluated by the recognition sensitivities, 
the overall recognition accuracy and the neurons needed. The recognition sensitivity to a particular pattern 
is defined as the ratio of beats of that pattern that are correctly recognized cto the total beats of that 
pattern, The overall accuracy is defined as the ratio of the total number of beats recognized correctly to 
the total number of beats in the test phase. Performance of the online BP network implemented with the 


g.BSAnalyzer is shown in table 4.4. 
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Figure 4.10: Jumper Function 


Network | N1 N2 Al A2 PVC1 | PVC2 Ll L2 Overall Accuracy 
BP [20] | 98% | 94% | 91.3% | 92% | 95.3% 96% 98.7% 100% 92.5% 
BP 90% | 93% | 94.5% | 88% | 90.71% 92% 89.6% | 95.79% 91.7% 
Table 4.4: Performance of the online network. 


4.4 Mobile Platform Implementation 


The algorithm consists of four major steps: 


e QRS detection 

e Template formation and adaptation 
e Feature extraction 

e Beat classification 


Arrhythmia Detector is an Android application for detection of abnormalities with the processing of 


ECG data. Intended to be used with wearable eHealth devices to gather empirical ECG data and also be 


connected to databases. 


4.4.1 Programming 


The software framework was implemented in JavaTM using the Android SDK 2.3.3 (Google Inc.). Android 


was used because of its open nature, widespread use and the portability of the code. Additionally it al- 
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Signal properties: Number of beats: Wave parameters: 


Sampling frequency: 360 Hz Normal: 2272 Mean QRS duration: 75 ms 
Number channels: 2 Complete bundle branche block: 1 Mean P duration: 144 ms 
Signal length Oh 30min 6s Incomplete bundle branche block: O Mean ST level: -0.05 m 
Recording date: Ventricular ectopic: o Mean PO interval: 197 ms 
Recording time: Supraventricular ectopic: 0 Mean QT interval: 384 ms 
Mean heart rate: NaN bpm Pause > 2.5 s: D Corrected OT interrablall ms 
Correlation class 1 Correlation class 3 Correlation class 5 
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Figure 4.11: Correlation Classes 


lowed simple integration of external ECG sensors via Bluetooth. The software framework consisted of three 


components: 


e A Data Delivery Service that provided data streaming from a Bluetooth connection or pre-recorded 


database files 
e A Signal Processing Service that implemented the above mentioned algorithms, and 


e A Graphical User Interface (GUI) that displayed the results. 


Data Delivery Service 


The Data Delivery Service allowed provisioning of ECG data to the actual processing algorithms in real- 
time. Implementation was realized as an Android background service. Using a service allowed to manage 


the Bluetooth connection to the sensor node independently of any foreground process, allowing the user to 
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Figure 4.12: Resting ECG 
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Figure 4.13: Driving ECG 


switch between applications without loosing the connection. Furthermore the service provided a consistent 
interface to the processing application for live and pre-recorded data. Data were either provided in live 
mode or database mode. In live mode, the service allowed acquisition of ECG lead II data via a Bluetooth 
connection to a computer privided of a database node. In database mode, the service supported reading of 
pre-recorded ECG lead II data from database files. In this case real-time sampling was simulated by the 


service. If available, the service additionally supported annotation characters for evaluation. 


Signal Processing Service 


The signal data provided by the Data Delivery Service were passed to the Signal Processing Service, which 


implemented the processing steps. To realize processing in real-time, the implementation was optimized with 
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Figure 4.14: Online Arrhtyhmia Classification 


respect to computational requirements and memory footprint. Circular buffers were used for all buffering 
operations during signal processing to avoid overhead and function calls were kept at a minimum. Al digital 


filters were automatically adapted to the sampling interval of the incoming data. 


Graphical User Interface 


A GUI was implemented, which allowed starting and stopping the processing and visualization of the results 
(Figure 4.15). To allow efficient plotting of signal data with high sampling rates, a specialized plotting 
component was implemented. Line-plots representing the raw ECG signal, extracted QRS complexes and 
variation of heart rate were displayed in the bottom area. Every time a QRS complex was detected and 
classified, it was marked either in green (normal beat) or red (abnormal beat). Additionally the current 
value of different features like heart rate, R-R interval in ms and number of recognized QRS complexes was 


displayed in the upper area. Upon closing an overview of the classification results was displayed. 


4.4.2 Evaluation 


For evaluation, the MIT-BIH Arrhythmia [13] databases were used. Using a PC, the records of both databases 
were downloaded and converted to a format readable by the application. The data were then uploaded to the 
memory card a mobile phone. The mobile phones used were: HTCTM Wildfire S A510e. For the evaluation, 
the currently valid annotation character supplied with the MIT-BIH records (if available) was compared with 
the result of the classification of every detected beat. The results were stored by the application to provide 
statistics about matches and mismatches. After processing of all MIT-BIH records, the results were evaluated 
on a PC. For the presentation of the results, the normal/abnormal beat classifications were evaluated. The 


measures to describe the performance of the classification were: 


67 


© 
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Figure 4.15: Graphical User Interface 


e TRUE NEGATIVE: correctly classified as normal 

e TRUE POSITIVE: correctly classified as abnormal 

e FALSE NEGATIVE: incorrectly classified as normal 
e FALSE POSITIVE: incorrectly classified as abnormal 


e MEDIAN FN: median of false negatives over all records 


Table 4.5 presents the results for real-time classification of the MIT-BIH Arrhythmia. 165, 898 in 72 records) 
were processed, and 0.42% were not recognized. Overall sensitivity for abnormal beat detection was 89.5% 
with a specificity of 80.6%. The results of the application were identical on the employed mobile phone and 


simulated ones working over Android Studio. 
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Parameter Measurement 
detected Beats 99.59% 
True Positive 11224 
True Negative 65855 
False Positive 10987 
False Negative 1680 

Median FN 3 


Table 4.5: Performance of the mobile implementation 
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Chapter 5 


Conclusions 


During the development of this work the following conclusions were found: 


e Asa result of an experimental phase in which QRS complexes and the main ECG characteristics were 
found, the best combination of features was selected. Compared to previous works for arrhythmia de- 
tection used around 21 [20] features creating a very complex wall of multiplications, this study achieved 
to have similar (91,7 % instead of 92,5%) overall accuracy when classifying the samples with a selec- 
tion of only 8 different values extracted from the QRS complexes. Similarly, the performance of 
three different neural networks was compared in order to discover that Back Propagation (BP) Neural 
Network was only 4% over Learning Vector Quantization (LVQ) Neural Network and above 3% of Self 
Organizing Maps (SOM) Neural Network. 


e The usage of Neural Network toolbox in Matlab lead to the conclusion that there is not such a big 
difference between these three types of training when the overall performance is taking into account. 
However, when talking about complexity, Back Propagation has the best balance between the algo- 
rithm complication and accuracy. The fact that BP training method is easier to implement also lead 


the study to conclude that this is the best choice for a mobile platform implementation. 


e After the work of training the BP network, all the values for the weights were found and posteriorly im- 
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plemented for classifying the data. 91,7% of times the trained network classified the samples correctly. 
Additionally, new samples were collected using gTech sensors and amplifiers for the instrumentation, 


to obtain an overall accuracy of 96,65% in this classifications. 


e Implementation in an Android mobile platform was achieved as a project of different stages: Con- 
nectivity, User Interface and also Performance of the network were tested. The experimental results 
showed that the complexity of the network can be maintained low without sacrificing the accuracy. 
Supplementary, other types of alert were implemented in the network since the connection with the 


car allows to create sound alarms and controls for the driver in real time. 


5.1 Future Work 


Posterior to this study, we propose future works related to other subsystems of the arrhythmia detector and 


also in the same line of software implementation: 


e Actual implementation of the system consumes lots of resources from the mobile device including en- 
ergy source. One possible solution is to develop a embedded system which implements the window 


extraction and the arrhythmia classification. 


e Possible linking with a database with physicians on the other side can be implemented for a deeper 


analysis and medical decisions. 


e Methodology for selection of features and training the network can be implemented for other types 
of biopotentials such as Electroencephalogram (EEG) or Electromiogram (EMG) signals to reach a 


better and more complete assistance environment. 


e Instrumentation has a very interesting opportunity of research in order to connect the contact-less 


sensors to the Bluetooth modules of the mobiles platforms. 
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e After the classification is done, the car can be also controlled by means of its actuators and warning 


signals in order to prevent possible fatal situations. 
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Appendix A 


ECG Filtering Code 


YH load ecg: MITDB 

[tm, x]J=rdsamp('mitdb/115' ,[] ,1000); 

plot (tm,x(:,1)) 

Fs=1/tm; 

%#: plot noisy signal 

figure 

subplot (911), plot(x), set(gca, 'YLim', [-1 1], 
title ('noisy') 


7H moving average filter 


window = 15; 
h = ones(window,1)/ window; 
y = filter(h, 1, x); 


subplot (914), plot(y), set(gca, 'YLim', [-1 1], 
title ('moving average') 


XÆ moving weighted window 


window = 7; 
h = gausswin(2* window+1)./window; 
y = zeros(size(x)); 


for i=l:length (x) 


for j=-window: window; 
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'xtick' ,[]) 


'xtick' ,[]) 


if j>-i && j<(length(x)-i+1) 
Ayi) = y(i) + x(i4j) + (1-(j/window)*2) / window; 
y(i) = y(i) + x(i+j) * h(j+window+1); 
end 
end 
end 
subplot (915), plot( y ), set(gca, 'YLim', [-1 1], 'xtick' ,[]) 
title ('weighted window') 
YH gaussian 


window = 7; 


h = normpdf( - window: window, 0, fix((2*window+1)/6) ); 

y = filter(h, 1, x); 

subplot(916), plot( y ), set(gca, 'YLim', [-1 1], 'xtick' ,[]) 
title ('gaussian') 

YH filter 

order = 15; 


h = firl(order, 0.1, rectwin(order+1)); 


II 


y filter(h, 1, x); 

subplot (918), plot( y ), set(gca, 'YLim', [-1 1], 'xtick' ,[]) 
title('firl') 

% lowpass Butterworth filter 

fNorm = 25 / (Fs/2); X% normalized cutoff frequency 
[b,a] = butter(10, fNorm, 'low'); %# 10th order filter 

y = filtfilt (b, a, x); 

subplot (919), plot(y), set(gca, 'YLim', [-1 1]) 


title ('butterworth') 
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Appendix B 


R-peak location 


function [pks,locs] = findRPeaks(X, varargin) 


narginchk (1,11); 


X,Ph,Pd,Th,Np, Str ,infldx] = parse_inputs(X, varargin {:}); 

pks ,locs] = getPeaksAboveMinPeakHeight (X, Ph); 

pks , locs] = removePeaksBelowThreshold (X, pks , locs ,Th, infIdx ); 

pks , locs] = removePeaksSeparatedByLessThanMinPeakDistance(pks,locs ,Pd); 
pks , locs] = orderPeaks (pks, locs , Str); 

pks , locs] = keepAtMostNpPeaks (pks , locs ,Np); 


function [X,Ph,Pd,Th,Np, Str ,infIdx] = parse_inputs(X, varargin) 
% Validate input signal 
validateattributes(X,f 'numeric'},{'nonempty','real','vector'},... 


'findpeaks','X'); 


try 
% Check the input data type. Single precision is not supported. 
chkinputdatatype(X); 

catch ME 
throwAsCaller (ME); 

end 
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error (message ( ' signal: findpeaks:emptyDataSet')); 
end 
Vitfunction dspopts.findpeaks 
hopts = uddpvparse('dspopts.findpeaks',varargin {:}); 
Ph = hopts.MinPeakHeight ; 
Pd = hopts.MinPeakDistance; 
Th = hopts. Threshold; 
Np = hopts.NPeaks; 


Str = hopts.SortStr; 


if ~isempty(Pd) && (~isnumeric(Pd) || ~isscalar(Pd) ||any(rem(Pd,1)) || (Pd < 1)) 


error (message ( ' signal: findpeaks:invalidMinPeakDistance', 


end 
% Set default values for MinPeakDistance and NPeaks 


if (isempty(Pd)), Pd 


II 
= 
o 
B 
a 


if (isempty(Np)), Np = M; end 
if (Pd >= M) 
error (message('signal: findpeaks: largeMinPeakDistance', 


end 


'MinPeakDistance')); 


'MinPeakDistance', 'MinPeakDistance', 


% Replace Inf by realmax because the diff of two Infs is not a number 


infldx = isinf(X); 
if any(infldx), 
X(infldx) = sign (X(infIdx ))*realmax; 
end 
infldx = infldx € X>0; % Keep only track of +Inf 
function [pks,locs] = getPeaksAboveMinPeakHeight (X, Ph) 
pks = []; 
locs = []; 
if all(isnan(X)), 
return, 
end 
Indx = find(X > Ph); 
if (isempty (Indx)) 
warning (message('signal: findpeaks:largeMinPeakHeight', 
return 


end 
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'MinPeakHeight' , 


'MinPeakHeight' )); 


I 


trend = sign(diff(X)); 
idx = find(trend==0); % Find flats 
N = length (trend); 
for i=length(idx):-1:1, 
% Back-propagate trend for flats 
if trend(min(idx(i)+1,N))>=0, 
trend(idx(i)) = 1; 
else 
trend(idx(i)) = -1; % Flat peak 
end 
end 
function [pks,locs] = removePeaksBelowThreshold (X, pks , locs ,Th,infIdx) 
idelete = []; 
for i = 1:length (pks), 
delta = min(pks(i)-X(locs(i)-1),pks(i)-X(locs(i)+1)); 
if delta<Th, 
idelete = [idelete i]; %fok<AGROW> 
end 
end 
if ~isempty (idelete), 


locs (idelete) = []; 


end 
X(infIdx) = Inf; % Restore +Inf 
locs = union(locs ,find(infldx)); % Make sure we find peaks like [realmax Inf realmax] 


pks =X(locs); 
function [pks,locs] = removePeaksSeparatedByLessThanMinPeakDistance(pks ,locs ,Pd) 
% Start with the larger peaks to make sure we don't accidentally keep a 
% small peak and remove a large peak in its neighborhood. 
if isempty(pks) || Pd==1, 
return 
end 
% Order peaks from large to small 
[pks, idx] = sort(pks,'descend'); 
locs = locs (idx); 


idelete = ones(size(locs)) <0; 
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for i = 1:length (locs), 
if ~idelete(i), 
% If the peak is not in the neighborhood of a larger peak, 
% secondary peaks to eliminate. 
idelete = idelete | (locs>=locs(i)-Pd)&(locs<=locs(i)+Pd); 
idelete(i) = 0; % Keep current peak 
end 
end 
pks(idelete) = []; 
locs(idelete) = []; 
function [pks,locs] = orderPeaks(pks,locs ,Str) 
if isempty(pks), return; end 
if stremp(Str,'none') 
[locs ,idx] = sort(locs); 
pks = pks(idx); 
else 
[pks,s] = sort(pks,Str); 
locs = locs(s); 
end 
function [pks,locs] = keepAtMostNpPeaks(pks,locs ,Np) 
if length (pks)>Np, 
locs = locs(1:Np); 
pks = pks(1:Np); 


end 
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Appendix C 


ECG Filtering Code 


function [qrs-amp-raw ,qrs-i-raw , delay]|=windowFinding (ecg, fs , gr) 
if ~isvector (ecg) 
error('ecg must be a row or column vector'); 
end 
if nargin < 3 


gr= i1; % on default the function always plots 
ecg = ecg(:); % vectorize 


%% Initialize 


qrs_c =|]; %amplitude of R 


qrs-i =[]; “index 

SIG_LEV = 0; 

nois_c =|]; 

nois_i =Í]; 

delay = 0; 

skip = 0; % becomes one when a T wave is detected 
not_nois = 0; % it is not noise when not_nois = 1 
selected_RR =[]; % Selected RR intervals 


m-_selected-RR = 0; 
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mean_RR = 0; 

qrs_i_raw =Í]; 
qrs_amp-_raw =|]; 
ser_back = 0; 

test_m = 0; 

SIGL_buf = []; 
NOISL_buf = ; 
THRS_buf = []; 
SIGL_bufl = |]; 
NOISL_bufl = []; 


THRS-_bufl = []; 
Y% Plot differently based on filtering settings 
if gr 
if fs == 200 
figure, ax(1)=subplot (321); plot(ecg);axis tight; title ('Raw ECG Signal'); 
else 
figure, ax(1)=subplot(3,2,[1 2]); plot(ecg);axis tight; title ('Raw ECG Signal'); 
end 
end 
%% Noise cancelation (Filtering) % Filters (Filter in between 5-15 Hz) 
if fs = 200 
Y Low Pass Filter H(z) = ((1 - z*(-6))"2)/(1 - 2*(-1))°2 
b= [100000 -200000 1]; 
a= [1 -2 1]; 
h-1 = filter(b,a,[1 zeros(1,12)]); 
ecg-l = conv (ecg ,h-_l); 
ecg-l = ecg-1/ max( abs(ecg-1)); 
delay = 6; %based on the paper 
if gr 
ax(2)=subplot(322);plot(ecg-1);axis tight;title('Low pass filtered'); 
end 
7% High Pass filter H(z) = (-14+32z" (-16)+z” (-32))/(1+z”*(-1)) 
b= [-1 000000000000000 32 -3200000000000000 1]; 
a = |1 -1]; 


h-h = filter (b,a,[1 zeros(1,32)]); 
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ecg-h = conv (ecg-1 ,h_h); 

ecg-h = ecg-h/ max( abs(ecg-h)); 

delay = delay + 16; % 16 samples for highpass filtering 

if gr 

ax(3)=subplot(323);plot(ecg-h);axis tight;title('High Pass Filtered'); 
end 

else 

Y% bandpass filter for Noise cancelation of other sampling frequencies(Filtering) 
f1=5; %cuttoff low frequency to get rid of baseline wander 

f2=15; %cuttoff frequency to discard high frequency noise 

Wn=[fl f2]*2/fs; % cutt off based on fs 

N = 3; % order of 3 less processing 

[a,b] = butter(N,Wn); %bandpass filtering 

ecg_h = filtfilt (a,b,ecg); 

ecg-h = ecg_h/ max( abs(ecg-h)); 

if gr 

ax(3)=subplot (323); plot(ecg_h);axis tight;title('Band Pass Filtered'); 
end 

end 

X% derivative filter H(z) = (1/8T)(-z*(-2) - 2z*(-1) + 2z + 2°(2)) 
hd = [-1 -2 0 2 1]*(1/8);%1/8* fs 

ecg_d = conv (ecg_h ,h_d); 

ecg_d = ecg_d/max(ecg_d); 

delay = delay + 2; % delay of derivative filter 2 samples 

if gr 

ax(4)=subplot (324); plot(ecg_d);axis tight;title('Filtered with the derivative filter'); 
end 

X% Squaring nonlinearly enhance the dominant peaks 

ecg_s = ecg_d.°2; 

if gr 

ax(5)=subplot (325); plot(ecg_s);axis tight;title('Squared'); 


end 


84 


Appendix D 


Training Algorithm 


close all, clear all, cle, format compact 
% generate data 


[X, Xtrain , Ytrain, fig] = data_generator (); 


% no hidden layers 


net = feedforwardnet ([]); 


Design 


% % one hidden layer with linear transfer functions 


% net = feedforwardnet ([10]); 

% net.layers{1l}.transferFcn = 'purelin'; 

% set early stopping parameters 
net.divideParam.trainRatio = 1.0; % training 
net .divideParam.valRatio = 0.0; % validation 


net.divideParam.testRatio = 0.0; % test set 


% train a neural network 


net.trainParam.epochs = 200; 
net = train (net , Xtrain , Ytrain); 
Ook RR E NO 


% view net 


view (net) 


set [%] 
set [%] 
[%] 


% simulate a network over complete input range 
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Y = net(X); 

% plot network response 

figure(fig) 

plot(X,Y,'color',[1 .4 0]) 

legend('original function','available data','Linear regression','location','northwest') 
% generate data 


[X, Xtrain , Ytrain , fig] = data_generator (); 


% choose a spread constant 
spread = .4; 
% create a neural network 


net = newrbe(Xtrain , Ytrain ,spread ); 


% view net 

view (net) 

% simulate a network over complete input range 

Y = net(X); 

% plot network response 

figure(fig) 

plot(X,Y,'r') 

legend ('original function','available data','Exact RBFN','location' ,'northwest') 
% generate data 


[X, Xtrain , Ytrain,fig] = data_generator (); 


% choose a spread constant 
spread = .12; 
% create a neural network 


net = newgrnn(Xtrain , Ytrain ,spread ); 


% view net 

view (net) 

% simulate a network over complete input range 
Y = net(X); 

% plot network response 


figure(fig) 
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plot(X,Y,'r') 
legend ('original function','available data','RBFN','location' ,'northwest') 
% generate data 


[X, Xtrain , Ytrain, fig] = data_generator (); 


% choose a spread constant 
spread = .2; 

% choose max number of neurons 
K = 20; 

% performance goal (SSE) 

goal = 0; 


% number of neurons to add between displays 


% create a neural network 


net = newrb(Xtrain , Ytrain , goal ,spread ,K, Ki); 


view (net) 

% simulate a network over complete input range 

Y = net (X); 

% plot network response 

figure(fig) 

plot (X,Y, Y) 

% Show RBFN centers 

c = net.iw{1}; 

plot(c,zeros(size(c)),'rs') 

legend('original function','available data','RBFN','centers','location','northwest') 
nono trainbr --------------- 

% Retrain a RBFN using Bayesian regularization backpropagation 
net.trainFen='trainbr'; 

net .trainParam.epochs = 100; 

% perform Levenberg-Marquardt training with Bayesian regularization 


net = train (net , Xtrain , Ytrain); 


% simulate a network over complete input range 


Y = net (X); 
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% plot network response 
figure(fig) 

plot (X,Y, 'm') 

% Show RBFN centers 

c = net.iw{1}; 


plot(c,ones(size(c)),'ms') 


legend('original function','available data','RBFN','centers','RBFN + trainbr' 


% generate data 


[X, Xtrain , Ytrain , fig] = data_generator (); 


% create a neural network 

net = feedforwardnet ([12 6]); 

% set early stopping parameters 
net.divideParam.trainRatio = 1.0; % training set [%] 
net.divideParam.valRatio = 0.0; % validation set [%] 
net .divideParam.testRatio = 0.0; % test set [%] 


% train a neural network 


net.trainParam.epochs = 200; 
net = train (net , Xtrain , Ytrain); 
A RN A este ete 


view (net) 

% simulate a network over complete input range 
Y = net(X); 

% plot network response 

figure (fig) 

plot(X,Y,'color',[1 .4 0]) 


legend('original function','available data','MLP','location 
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